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ABSTRACT 

Context. Theoretical simulations and observations at different angular resolutions have shown that magnetic fields have a central role 
in massive star formation. Like in low-mass star formation, the magnetic field in massive young stellar objects can either be oriented 
along the outflow axis or randomly. 

Aims. Measuring the magnetic field at milliarcsecond resolution (10-100 au) around a substantial number of massive young stellar 
objects permits determining with a high statistical significance whether the direction of the magnetic field is correlated with the 
orientation of the outflow axis or not. 

Methods. In late 2012, we started a large VLBI campaign with the European VLBI Network to measure the linearly and circularly 
polarized emission of 6.7 GHz CH3OH masers around a sample of massive star-forming regions. This paper focuses on the first seven 
observed sources, G24.78-l-0.08, G25.65-HL05, G29.86-0.04, G35.03-l-0.35, G37.43-I-L5L G174.20-0.08, and G213.70-12.6. For all 
these sources, molecular outflows have been detected in the past. 

Results. We detected a total of 176 CH3OH masing cloudlets toward the seven massive star-forming regions, 19% of which show 
linearly polarized emission. The CH3OH masers around the massive young stellar object MMl in G174.20-0.08 show neither linearly 
nor circularly polarized emission. The linear polarization vectors are well ordered in all the other massive young stellar objects. We 
measured significant Zeeman splitting toward both A1 and A2 in G24.78-l-0.08, and toward G29.86-0.04 and G213.70-12.6. 
Conclusions. By considering all the 19 massive young stellar objects reported in the literature for which both the orientation of the 
magnetic field at milliarcsecond resolution and the orientation of outflow axes are known, we find evidence that the magnetic field 
(on scales 10-100 au) is preferentially oriented along the outflow axes. 

Key words. Stars: formation - masers: methanol - polarization - magnetic fields - ISM: individual: G24.78-l-0.08, G25.65-1-1.05, 
G29.86-0.04, G35.03 -h0.35, G37.43-hL51, G174.20-0.08, G213.70-12.6. 


1. Introduction 

core accretion model describes the formation of high-mass 
i as a scaled-up version of the formation process of low- 
5 stars (e.g., McKee & Tan 2003). Specifically, it is proposed 
that massive stars form through gravitational collapse, which 
involves disk-assisted accretion to overcome radiation pressure 
and matter-ejection perpendicular to the disk to redistribute the 
angular momentum (e.g., McKee & Tan 2003). Nevertheless, 
only when the magnetic field has been taken into consideration, 
the theoretical simulations begin to faithfully reproduce the ob¬ 
servations (e.g., Peters et al. 2011; Seifried et al. 2012; Myers 
et al. 2013). This suggests that magnetic fields might play a role 
in massive star formation just as importantly as in the formation 
of low-mass stars. In low-mass star formation the magnetic field 
is thought to slow the collapse, to transfer the angular momen¬ 
tum, and to power the outflow (e.g., McKee & Ostriker 2007). 
However, in studies on low-mass star formation it is still an 


open debate whether the magnetic field aligns with the molec¬ 
ular outflows. Recently, two independent polarization surveys 
of low-mass protostellar cores, which were carried out at dif¬ 
ferent spatial resolutions, showed two contrasting results. Hull 
et al. (2013) found on scales of a few 100 to 1000 au no cor¬ 
relation between magnetic field orientation and outflow axis in 
low-mass young stellar objects (YSOs), while Chapman et al. 
(2013) found a good alignment on larger scales (> 2 x 10^ au). 
Similar conflicting results were also found toward massive YSOs 
(Surcis et al. 2012, 2013; hereafter Paper I and Paper II, re¬ 
spectively; Zhang et al. 2014). Surcis and collaborators started 
a VLBI-observations campaign of 6.7 GHz CH 3 OH masers to¬ 
ward massive star-forming regions (SFRs) to determine if there 
exists any correlation between the orientations of the magnetic 
field and of the outflow on very small scales (tens of au). Based 
on the results of nine sources, they found evidence that on scales 
of 10-100 au the magnetic field around massive YSOs is prefer¬ 
entially oriented along the outflow (Paper II). On the other hand. 
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based on a larger sample (21 sources), Zhang et al. (2014) re¬ 
ported that at arcsecond resolution (thousands of au) the outflow 
axis appears to be randomly oriented with respect to the mag¬ 
netic field in the core. These contrasting results might be due to 
the different resolutions; indeed, Zhang et al. (2014) postulated 
that angular momentum and dynamic interactions, possibly due 
to close binary or multiple systems, dominate magnetic fields at 
scales of about 10 ^ au. 

Before drawing any conclusion, it is important to improve 
the statistics by enlarging the number of massive SFRs toward 
which the orientation of the magnetic field at milliarcsecond 
(mas) resolution has been measured. Therefore, we have se¬ 
lected a flux-limited sample of 37 massive SFRs with declination 
> -9° and a total CH 3 OH maser single-dish flux greater than 
50 Jy from the 6.7 GHz CH 3 OH maser catalog of Pestalozzi et 
al. (2005). To increase the likelihood of detecting circularly po¬ 
larized CH 3 OH maser emission (< 1%) and thus allow the de¬ 
termination of the magnetic field strength, we have excluded the 
six regions hosting CH 3 OH maser that in recent single-dish ob¬ 
servations showed a total flux below 20 Jy (Vlemmings et al. 
2011). The total number of massive SFRs of the flux-limited 
sample is thus 31. The polarimetric 6.7 GHz CH 3 OH maser 
observations, and the subsequent measurement of the magnetic 
field orientation, of twelve of these SFRs had already been pub¬ 
lished in the recent past (Vlemmings et al. 2010; Surcis et al. 
2009, 2011a, 2014; Papers I and II). Therefore, 19 massive SFRs 
remain to be observed. We were given European VLBI Network^ 
(EVN) time to observe all of them at 6.7 GHz in several sessions 
between November 2012 and June 2015 (see Sect. 3). Here, we 
present the results of the first seven observed sources. The re¬ 
sults of the remaining twelve sources will be published in future 
papers of the present series as soon as they are observed and the 
data are fully analyzed. Eollowing the organization of Paper I 
and II, the sources are briefly introduced in Sects. 2.1-2.7, while 
the observations and our analysis are described in Sect. 3. The 
results, which are presented in Sect. 4, are discussed in Sect. 5, 
where we update our previous statistics. 

2. Massive star-forming regions 

2.1. G24.78+0.08 

G24.78-I-0.08 is one of the most studied massive SER (e.g., 
Codella et al. 1997; Cesaroni et al. 2003; Beltran et al. 2006, 
2011). The region is located at a kinematic distance of 7.7 kpc 
(Codella et al. 1997) and contains four centers of star formation, 
named from A to D (Euruya et al. 2002). Clump A is composed 
of two subclumps (Al and A2) that had been resolved in five dis¬ 
tinct cores by Beltran et al. (2011). These cores (named Al, Alb, 
Ale, A2, and A2b) are aligned in a southeast-northwest direc¬ 
tion coincident with the CO-outflow (position angle = 

-40°) that is associated with core A2 (Mai = 22 Mq; Beltran 
et al. 2011). Codella et al. (2013) confirmed the source embed¬ 
ded in A2 as the driving source of the outflow by imaging the 
SiO emission (PA^^°g^^ ss -45°) at arcsecond resolution with 
the SMA. Core Al (Mai - 16 Mq; Beltran et al. 2011) is as¬ 
sociated with a hypercompact (HC) H ii region (Galvan-Madrid 
et al. 2008). Both Al and A2 are embedded in toroids that ro¬ 
tate clockwise with position angle PAai = -1-50° and PAa 2 = 
-1-40° (Beltran et al. 2004, 2005, 2011), that is, the toroids are al¬ 
most perpendicular to the axis of the CO-outflow. Moscadelli 

' The European VLBI Network is a joint facility of European, Chinese, 
South African and other radio astronomy institutes funded by their na¬ 
tional research councils. 


et al. (2007) detected 6.7 GHz CH 3 OH masers around both 
Al and around A2. The CH 3 OH masers show velocity distri¬ 
butions consistent with the velocities of the toroidal structures 
(^toroid = 109.0-112.2kms-' and ^ = 109.6-lll.6kms-'; 
Beltran et al. 2011). Moreover, Moscadelli et al. (2007) also de¬ 
tected 22 GHz H 2 O masers around Al. These masers trace a fast 
(~40 km s ') expanding shell surrounding the HC H ii region. 

Using observations made with the Effelsberg 100 m tele¬ 
scope, Vlemmings et al. (2011) measured a Zeeman splitting of 
the CH 3 OH maser emission of AVz = (h- 0.50 + 0.08) m s '. 


2.2. G25.65+1.05 

The massive SER G25.65-i-l.05 (also known as IRAS 18316- 
0602 and RAEGL7009S) is located at a kinematic distance 
of 3.17 kpc (Molinari et al. 1996). The region is associated 
with a weak and irregular compact radio source that was ini¬ 
tially classified as an ultracompact (UC) H ii region (Kurtz et al. 
1994; Walsh et al. 1998). The radio source spatially coincides 
with an unresolved infrared source (Zavagno et al. 2002; Varri- 
catt et al. 2010) and with submillimeter emissions at 350 pm, 
450 pm, and 850 pm (Hunter et al. 2000; Walsh et al. 2003). 
A bipolar CO-outflow (PA^j^'j'^^ ss -65°) centered on the ra¬ 
dio source was first detected by Shepherd & Churchwell (1996). 
Recently, Sanchez-Monge et al. (2013) mapped the outflow us¬ 
ing a more reliable jet tracer, SiO emission. They detected both 
the red- (h- 45.8 kms-'< < - 1 - 88 .1 kms-') and blue- 

shifted (-1-5.9 kms-'< ^ -1-39.5 kms"') lobes of the 

jet or outflow (PA®^^^^ = -15°). Eurthermore, four 6.7 GHz 
CH 3 OH masers were detected near the continuum peak of the ra¬ 
dio source; they are linearly distributed southward (Walsh et al. 
1998). The CH 3 OH maser velocities suggest an association with 
the radio source, possibly with a disk and not with the bipolar 
outflow (Zavagno et al. 2002). 

Einally, Vallee & Bastien (2000) mapped the magnetic field 
toward the radio source at 760 pm, finding an orientation of 
the magnetic field of = - 1 - 8 ° + 16° (scale of lO"* au). A 

Zeeman splitting of the 6.7 GHz CH 3 OH maser emission of 
AVz = (-1-0.46 + 0.05) ms-' was measured with the Effelsberg 
100 m telescope (Vlemmings et al. 2011). 


2.3. G29.86-0.04 

G29.86-0.04 is at a kinematic distance of 7.4 kpc and has a 
velocity = H-101.85 kms-' (de Villiers et al. 2014). 

Caswell et'al. (1993, 1995) detected 12 GHz and 6.7 GHz 
CH 3 OH masers toward the region. The 6.7 GHz CH 3 OH masers 
show an arched distribution (150 mas x 340 mas) accompa¬ 
nied by a clear velocity gradient at mas resolution (Eujisawa 
et al. 2014). The CH 3 OH masers are associated with one of 
the two cores that were detected toward the region (Hill et al. 
2005, 2006). No 22 GHz H 2 O masers have been detected (Breen 
& Ellingsen 2011). A bipolar CO-outflow is associated with 
the 6.7 GHz CH 3 OH masers (de Villiers et al. 2014). While 
the redshifted lobe (-1-104 kms"'< < H-110 kms"') of 

the outflow is oriented almost south-north on the plane of the 
sky ~ - 1 - 6 °), the blueshifted lobe (-1-90 kms"'< 

^ -1-96.5 kms"') bends westwards passing from about 

6 ° to ^ +60° (de Villiers et al. 2014). 

The 6.7 GHz CH 3 OH maser Zeeman-splitting was measured 
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Table 1. Observational details. 


(1) 

Source 

(2) 

Program 

code 

(3) 

observation 

date 

(4) 

Calibrator 

(5) 

Polarization 

angle 

C) 

(6) 

Beam size 

(mas X mas) 

(7) 

Position 

Angle 

o 

(8) 

rms 

,_n^. 

(9) 

O's.-n.'’ 

( mJy s 

(10) (11) (12) (13) 

Estimated absolute position using FRMAP 
^2000 ^2000 A6“ 

(^ : ™ (mas) (mas) 

G24.78+0.08 

ES072 

30 May 2013 

J2202+4216 

-31 ±4 

10.4 X 4.0 

-36.14 

4 

8 

-(•18:36:12.563 

-07:12:10.787 

0.4 

3.7 

G25.65+1.05 

ES072 

31 May 2013 

J2202+4216 

-31 ±4 

11.9x3.5 

-39.94 

2 

35 

+ 18:34:20.900 

-05:59:42.098 

2.3 

18.3 

G29.86-0.04 

ES072 

01 June 2013 

J2202+4216 

-31 ±4 

9.0 X 3.6 

-40.71 

4 

6 

+ 18:45:59.572 

-02:45:01.573 

8.3 

172.5 

G35.03+0.35 

ES069 

04 Nov. 2012 

J2202+4216 

-30 ±2 

6.3 X 4.8 

-34.53 

4 

6 

+ 18:54:00.660 

+02:01:18.551 

7.2 

167.7 

G37.43+1.51 

ES072 

02 June 2013 

J2202+4216 

-31 ±4 

8.3 X 3.8 

-52.43 

3 

22 

+ 18:54:14.229 

+04:41:41.138 

7.0 

81.2 

G174.20-0.08 

ES069 

04 Nov. 2012 

J0555+3948 

-73 ±5 

7.7 X 4.3 

-28.78 

4 

5 

+05:30:48.020 

+33:47:54.611 

0.7 

1.0 

G213.70-12.6 

ES069 

03 Nov. 2012 

J0555+3948 

-73 ±5 

7.4 X 5.4 

-3.13 

4 

10 

+06:07:47.860 

-06:22:56.626 

2.1 

17.9 


Notes. Formal errors of the fringe rate mapping. Self-noise in the maser emission channels (e.g., Sault 2012). 


to be AVz = (-1-0.50 + 0.08) ms ' with the Effelsberg 100 m 
telescope (Vlemmings et al. 2011). 

2.4. G35.03+0.35 

The extended green object (EGO) G35.03 h- 0.35 hosts sev¬ 
eral massive YSOs at early evolutionary stages (Cyganowski 
et al. 2009; Paron et al. 2012). This massive SER (Visr = 
-1-51.5 kms^'; Paron et al. 2012) is located at a kinematic dis¬ 
tance of 3 . 43 j^Q 3 g kpc (Cyganowski et al. 2009). Pour of the 
five radio continuum sources that were detected toward the re¬ 
gion (CMl-5) are aligned with the bipolar morphology of the 
4.5 jum emission (PA”*-^^™ = +27°; Cyganowski et al. 2011). 
CMl, which is a well-known UC H ii region, and CM4 are asso¬ 
ciated with the southwestern lobe of the 4.5 //m emission, CM3 
is associated with the northeastern lobe, and CM2 is located be¬ 
tween the two lobes. The symmetric spacing of CM3 and CM4 
relative to CM2 might be the signature of knots in an ionized 
jet (Cyganowski et al. 2011). Purthermore, the radio spectral 
index of CM2 suggests that the radio source might either be a 
HC H II region or the product of an ionized wind that hits the 
surrounding gas (Cyganowski et al. 2011; Paron et al. 2012). 
Paron et al. (2012) detected a bipolar ^^CO-outflow at a resolu¬ 
tion of tens of arcseconds (beam size = 22 arcsec) that is coin¬ 
cident in position with the whole 4.5 fim emission. The axis of 
the ^^CO-outflow is oriented almost along the line of sight, with 
the redshifted lobe (-1-58 kms“'< < +66 kms“') south¬ 
east and the blueshifted lobe (+37 km s“^ < < +49 km s“') 

northwest of the axis. 

CH 3 OH, H 2 O, and OH masers were detected toward CM2 
(Porster & Caswell 1999; Argon et al. 2000; Cyganowski et 
al. 2009; Pandian et al. 2011). The 6.7 GHz CH 3 OH masers, 
which are all blueshifted, lie on the “waist” between the two 
lobes of the 4.5 //m emission and show a complex morphol¬ 
ogy at scales of 10 mas (Cyganowski et al. 2009; Pandian et 
al. 2011). Recently, Caswell et al. (2013) measured a persistent 
linearly polarized emission of the OH masers over several years. 
A large Zeeman splitting of AVz = (+1.22 + 0.23) ms^'of the 
6.7 GHz CH 3 OH maser emission was measured by Vlemmings 
et al. (2011) with the Effelsberg 100 m telescope. 

2.5. G37.43+1.51 

The massive SER G37.43+1.51 coincides with the IRAS source 
18517+0437 (Visr = +44.1 kms^'; Lopez-Sepulcre et al. 2010), 
and it is located at a parallax distance of 1.88);QQg kpc (Wu et 
al. 2014). Lopez-Sepulcre et al. (2010) detected a C^^O-outflow 


oriented north-south = -4°) with the redshifted lobe 

(+46 kms“'< < +50 kms^') and the blueshifted lobe 

(+38 kms“'< Vy™ < +42 kms^') located north and south, 
respectively. The C'^O-outflow is associated with 6.7 GHz 
CH 3 OH masers that have been detected in a linear distribution 
northwest-southeast with a clear velocity gradient (Schutte et al. 
1993; Pujisawa et al. 2014; Wu et al. 2014). Vlemmings (2008) 
measured a Zeeman splitting of the 6.7 GHz CH 3 OH maser of 
AVz = (+0.75+ 0.09) ms-'. 

2.6. G174.20-0.08 

In G174.20-0.08, which is better known as APGL5142, 
two centers of massive star formation were identified: 
IRAS 05274+3345 and IRAS 05274+3345-East (Hunter et al. 
1995; Torrelles et al. 1992). This massive SER is located 
at a kinematic distance of 1.8 kpc (Snell et al. 1988). 
IRAS 05274+3345-East (Vur = -1.0 kms-'; Zhang et al. 2007) 
hosts five 1.3 mm cores (MM-1 to MM-5) and three CO- 
outflows (Zhang et al. 2007). Outflow-C is associated with core 
MM-1, which powers 22 GHz H 2 O and 6.7 GHz CH 3 OH masers 
(Goddi et al. 2007, 2011). While the proper-motion measure¬ 
ments of the H 2 O masers trace the expansion of the colli¬ 
mated outflow-C (PA“\°^^ = -40°; Goddi et al. 2011), the 
CH 3 OH masers instead trace an infall of gas onto the central 
massive protostar (Goddi et al. 2011). Palau et al. (2011) found 
evidence of a possible disk perpendicular to outflow-C by ob¬ 
serving complex organic molecules. No Zeeman splitting of the 
CH 3 OH maser emission was measured toward APGL 5142 with 
the Effelsberg 100 m telescope (< 0.08%; Vlemmings 2008). 


2.7. G213.70-12.6 

The source G213.70-12.6 (D = 0.83 Kpc; Herbst & Racine 
1976) is better known under the name Monoceros R2 (here¬ 
after Mon R2) and is composed of several H ii regions and YSOs 
(e.g., Howard et al. 1994; Carpenter et al. 1997; Preibisch et al. 
2002). G213.70-12.6 hosts several infrared sources, the brightest 
of which is IRS 3 (L = 14000 Lq; Henning et al. 1992). IRS 3 is 
a compact cluster of massive YSOs (Preibisch et al. 2002) that 
is located at ~ 50" northeast of the center of a giant CO-outflow 
(P^omflow ~ -45°) that is powered by IRS 6 (Xu et al. 2006; 
beam size 15"). Very recently, a bipolar *^CO(2- 1) outflow 
(PA^^*^^'^^*' = +53°) has been detected toward IRS 3 (Dierickx 
et al. 2015; beam size ~ 0.5"). Its redshifted lobe (+21 kms-'< 
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Fig. 1. Total intensity (I, upper panel) and circularly polarized intensity (VJ lower panel) spectra for the CH 3 OH maser features G24.16, G24.23, 
G24.43, G24.52, G29.09, and G213.15 (see Tables A.l, A.3, and A.7). The thick red lines are the best-fit models of 1 and V emission obtained 
using the adapted FRTM code (see Sect. 3). The maser features were centered on zero velocity. 


i3cof2—n _i 

< -1-26 kms ) is located southwest of the source and 

its blueshifted lobe (-5 kms“'< < -i-2 kms^') lies to 

the northeast. One of the massive YSOs of IRS 3, called Star A 
(12 Mo< M < 15 Mq; Preibisch et al. 2002), is associated with 
6.7 GHz CH 3 OH masers that lie in a northeast-southwest linear 
distribution of 170 mas (PAp™Qjj =» 40°; Minier et al. 2000). 
Moreover, Star A is located along the axis of the '^CO(2 - 1) 
outflow (Fig. 5 of Dierickx et al. 2015). 

Curran & Chrysostomou (2007) measured using polarimet- 
ric observations at 850 pm a magnetic field strength of ~0.2 mG 
throughout G213.70-12.6. The polarization percentage around 
IRS 3 decreases below 1%, and the magnetic field at a resolu¬ 
tion of 6'.'18 changes its orientation from north-south to east- 
west (see Fig. 1 of Curran & Chrysostomou 2007). At a spatial 
resolution of 0".97, the polarization vectors of the 2.16 pm emis¬ 
sion form an elliptical pattern in a ~ 10" region around IRS 3 
(PAB-pattern = -40°; Yao et al. 1997). More recently, Simpson 
et al. (2013) measured the near-infrared (2 pm) polarimetry of 
IRS 3 at a spatial resolution of 0"2. They found that the frac¬ 
tional linear polarization and the orientation of the linear polar¬ 
ization vectors around Stars A and B are consistent with the mea¬ 
surements at larger scale (Yao et al. 1997; Curran & Chrysosto¬ 
mou 2007), and they measured a polarization PA of Star A of 

p^Star A ^ _ggo ^ 2°. 


3. Observations and analysis 

The first seven massive SFRs were observed at 6.7 GHz in full 
polarization spectral mode with eight of the EVN antennas (Ef- 
felsberg, Jodrell, Onsala, Medicina, Noto, Torun, Westerbork, 
and Yebes-40m) between November 2012 and June 2013, for 
a total observation time of 49 h. The bandwidth was 2 MHz, 
providing a velocity range of ~ 100 kms '. The data were cor¬ 
related with the EVN software correlator (SEXC; Keimpema et 
al. 2015) at the Joint Institute for VLBI in Europe (JIVE) us¬ 
ing 2048 channels and generating all four polarization combi¬ 
nations (RR, LL, RE, LR) with a spectral resolution of ~ 1 kHz 
(~0.05 kms '). All the observational details are reported in Ta¬ 
ble 1. We report in Cols. 1 to 3 the target source, the program 
code, and the date of the observations; in Cols. 4 and 5 we 
list the polarization calibrators with their polarization angles. 
Columns 6 to 8 list the restoring beam sizes, corresponding po¬ 
sition angles, and the thermal noise. In Col. 9 we also show the 
self-noise in the maser emission channels (see below for more 
details). Einally, Cols. 10 to 13 report the estimated absolute po¬ 
sition of the reference maser and the ERMAP uncertainties (see 
below for more details). 

The data were edited and calibrated using AIPS. The band¬ 
pass, delay, phase, and polarization calibration were performed 
on the calibrators listed in Table 1. Eringe-fitting and self¬ 
calibration were performed on the brightest maser feature of 
each star-forming region. The /, Q, U, and V cubes were imaged 
using the AIPS task IMAGR. The Q and U cubes were combined 
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to produce cubes of polarized intensity (POLI - + U^) and 

polarization angle (POLA - 112 x atan{UIQ)). We calibrated 
the linear polarization angles by comparing the linear polariza¬ 
tion angles of the polarization calibrators measured by us with 
the angles obtained by calibrating the POLCAL observations 
made by NRAO^. The NRAO POLCAL observing program was 
temporarily interrupted because of the JVLA commissioning. 
The last POLCAL observations were made in May/June 2012, 
therefore we were able to calibrate the polarization angles of the 
sources observed in 2012 by using the results from the last ob¬ 
serving run. The calibrator observed in 2013 was J2202H-4216, 
which shows a constant polarization angle between 2005^ and 
2012 of -31°+4°. To calibrate the polarization angles of the 
maser sources observed in 2013, we therefore assumed that the 
polarization angle of J2202-I-4216 has not changed significantly. 
We were thus able to estimate the polarization angles with a sys¬ 
temic error of no more than ~ 5° (see Col. 5 of Table 1). The 
formal errors on POLA are due to thermal noise. This error is 
given by ctpola - 0.5 {crplPOLL) x (180°/7r) (Wardle & Kron- 
berg 1974), where <Tp is the rms error of POLL 

Because the observations were not performed in phase- 
referencing mode, we estimated the absolute position of the 
brightest maser feature of each source through fringe rate map¬ 
ping by using the AIPS task FRMAP. The results and the for¬ 
mal errors of FRMAP are reported in Cols. 10 to 13 of Table 1. 
The absolute positional uncertainties are dominated by the phase 
fluctuations that we estimate to be on the order of no more than 
a few mas from our experience with other experiments and vary¬ 
ing the task parameters. 

We analyzed the polarimetric data following the procedure 
reported in Papers I and IT First, we identified the CH 3 OH maser 
features by using the process described in Surcis et al. (201 lb), 
and then we determined the mean linear polarization fraction 
(Pi) and the mean linear polarization angle (x) across the spec¬ 
trum of each CH 3 OH maser feature. Second, we made use of the 
adapted full radiative transfer method (FRTM) code for 6.7 GHz 
CH 3 OH masers (Vlemmings et al. 2010, Surcis et al. 2011a, 
Paper II) to model the total intensity and the linearly polarized 
spectrum of every maser feature for which we were able to de¬ 
tect linearly polarized emission. The output of this code pro¬ 
vides estimates of the emerging brightness temperature (PbAQ) 
and of the intrinsic thermal line width (AV)). Following Surcis 
et al. ( 2011 a), we restricted our analysis to values of AVi from 
0.5 kms“^ to 1.95 km s '. From T\,AQ. and Pi, we then deter¬ 
mined the angle between the propagation direction of the maser 
radiation and the magnetic field (ff). If 0 > 0crit = 55°, where 0crit 
is the Van Vleck angle, the magnetic field appears to be perpen¬ 
dicular to the linear polarization vectors; otherwise, it is parallel 
(Goldreich et al. 1973). To better determine the orientation of 
the magnetic field with respect to the linear polarization vec¬ 
tors, we followed the method introduced in Paper II that takes 
into consideration the errors associated with 0, that is, Ac¬ 
cording to this, the magnetic field is most likely perpendicular to 
the linear polarization vectors if - 55°| > \0^ - 55°|, where 
0"^ = otherwise, the magnetic field is assumed to be paral¬ 

lel. Of course, if 0^ and 0^ are either larger or smaller than 55° , 
the magnetic field is perpendicular or parallel to the linear polar¬ 
ization vectors, respectively. 

Note that if T^Afl > 2.6 x 10® K sr the 6.7 GHz 

CH 3 OH masers can be considered partially saturated and their 
AVi and PbAQ values are overestimated and underestimated, re- 


^ http://www.aoc.nrao.edu/~smyers/evlapolcal/polcal_master.html 
^ http://www.aoc.nrao.edu/~smyers/calibration/ 


spectively (Surcis et al. 2011a). However, we are confident that 
the orientation of their linear polarization vectors is not affected 
by their saturation state (Paper I), and consequently, they can be 
used for determining the orientation of the magnetic field in the 
region. 

Finally, to measure the Zeeman splitting (AVz), we included 
the best estimates of TbAQ and AVi in the FRTM code to pro¬ 
duce the I and V models used for fitting the total intensity and 
circularly polarized spectra of the corresponding CH 3 OH maser 
feature (Fig. 1). Because the circularly polarized emission of 
CH 3 OH masers is usually very weak (< 1%), we must take 
into consideration the self noise"' (ctj ^n.) produced by the masers 
(Col. 9 of Table 1; e.g., Sault 2012) when we measure the Zee- 
man splitting. Therefore, we consider real a detection of circu¬ 
larly polarized emission only when the detected V peak flux 
of a maser feature is both five times higher than the rms and 
three times larger than ctj . We know from the Zeeman effect 
theory that AVz is related to the magnetic field strength along 
the line of sight (By) through AVz = az ■ By. However, the 
Lande g-factors for the CH 3 OH molecule (including the 6.7 GHz 
maser transition) on which az depends are still unknown, and 
consequently, the magnetic field strength cannot yet be derived 
from our Zeeman-splitting measurements (e.g., Vlemmings et al. 
2011). 

4. Results 

In Tables A. 1-A.7 we list all the 6.7 GHz CH 3 OH maser features 
detected toward the seven massive SFRs observed with the EVN. 
The description of the maser distribution and the polarization 
results are reported for each source separately in Sects. 4.1^.7. 
In Figs. 2-8 we show the measured linear polarization vectors as 
black segments and the inferred orientation of the magnetic field, 
which is either parallel or perpendicular to the linear polarization 
vectors (see Sect. 3), in green in the bottom right corner of each 
panel. 

4.1. G24.78+0.08 

We detected 53 CH 3 OH maser features, named G24.01-G24.53 
in Table A.l, 33 toward core Al and 20 toward A2. In Fig. 2 
we show all the maser features associated with Al in the left 
panel and those associated with A2 in the right panel. The 
maser distributions around the two cores are identical to those 
observed previously by Moscadelli et al. (2007), even though 
we detected about 40 maser features more. The peak flux 
density range (Col. 5 in Table A.l) and the local standard of rest 
velocity (Visr; Col. 6 in Table A.l) range are similar to previous 
measurements. 

We detected linearly polarized emission from ten 
CH 3 OH maser features around Al (BJ^' = 0.8% - 3.5%) 
and from three maser features around A2 (P^^ = 1.0% - 1.3%). 
The adapted FRTM code was able to fit all of them but G24.38. 
The outputs of the code are reported in Cols. 10, 11, and 14 of 
Table A.L The twelve maser features for which we estimated 
TbAQ are unsaturated. Indeed, T},AQ. < 2.6 X 10’' K sr (or 
in logarithmic value < 9.4 log K sr). For the maser features 
G24.23 and G24.52, both associated with Al, we have that 

- 55°| < \0^ - 55°|, that is, the magnetic field is assumed to be 
parallel to their linear polarization vectors as described in Sect. 3. 

The self-noise is high when the power contributed by the astronom¬ 
ical maser is a significant portion of the total received power (Sault 
2012 ). 
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Fig. 2. View of the CH3OH maser features detected around the mm subcores A1 (left panel) and A2 (right panel) of G24.78+0.08. The reference 
position is a '2000 = +18''36“12®.563 and (^2000 = -07°12'10"787 (see Sect. 3). Triangle symbols identify CH3OH maser features scaled logarith¬ 
mically according to their peak flux density (Table A. 1). Maser LSR radial velocities are indicated by color (the assumed velocity of the region 

I2r''r\t 1 _m ■ . 

^isr - -I-111 km s“ , Furuya et al. 2002). A 1 Jy bearn^ symbol is plotted for illustration in both panels. The linear polarization vectors, 

scaled logarithmically according to polarization fraction Fj, are overplotted. In the bottom right corner of both panels the corresponding error- 
weighted orientation of the magnetic field ((®b), see Sect.5.1) is also reported, the two dashed segments indicate the uncertainty. The two arrows 
in the right panel indicate the direction, and not the actual position, of the red- and blueshifted lobes of the *^CO(l - 0) outflow associated with 
G24.78-t-0.08-A2 (PA^^™^ = -40°; Beltran et al. 2011). The dotted lines indicate the direction of the CH 3 CN toroids (Beltran et al. 2011). The 
circled plus and minus symbols indicate where the magnetic field points away from and where toward the observer, respectively. 


We also measured AVz for five CH 3 OH maser features (Col. 13 
of Table A.l), only one of which is associated with A2. The 
circular polarization fraction (Py) ranges from 0.3% to 0.7% and 
the Zeeman splitting is -9.7 ms“*< < h- 7.8 ms“^ around 

Al and AV^^ = (-4.0 + 0.8) ms“' around A2 (see Fig. 1). 

4.2. G25.65+1.05 

Imaging a 2" x 2" field-of-view centered on G25.02, we were 
able to detect a total of 23 6.7-GHz CH 3 OH maser features, 
named G25.01-G25.23 in Table A.2. The maser features can be 
divided into two groups (named here group A and group B) sep¬ 
arated from each other by about 400 mas (~1300 au; see Fig. 3). 
The two groups are located at the origin of the bipolar outflow. 
Comparing our detections with the four CH 3 OH maser spots de¬ 
tected by Walsh et al. (1998), which were linearly distributed 
southwards over 1", we note that only group A can be associ¬ 
ated with one of the previous maser spots, spot B (as named by 
Walsh et al. 1998). The other three CH 3 OH maser spots were 
not detected by us, and group B was not detected by Walsh et al. 
(1998). All the maser features of group B are blueshifted with 
respect to the systemic velocity = -1-42.41 kms“'; 

Sanchez-Monge et al. 2013). Even though the maser features of 
group A show both a linear distribution (PAchsOh = -49° + 7°) 
and red- and blueshifted velocities, no clear velocity gradient is 
observed. 

Five CH 3 OH maser features of group A show linearly polar¬ 
ized emission (Pi = 0.3% -1.3%), and according to the output of 
the adapted FRTM code all of them are unsaturated (see Col. 11 
of Table A.2). For the maser features G25.02 and G25.06 we de¬ 
termined that - 55°| < \6^ - 55°|, that is, the magnetic field 
is assumed to be parallel to their linear polarization vectors. For 



Right Ascension offset (mas) 

Fig. 3 . View of the CH 3 OH maser features detected around 
G25.65-I-1.05 (Table A.2). Same symbols as in Fig. 2. The assumed ve¬ 
locity of the YSO is - -1-42.41 km s“* (Sanchez-Monge et al. 

2013). The two arrows indicate the direction, and not the actual position, 
of the red- and blueshifted lobe of the bipolar outflow = -15°; 

Sanchez-Monge et al. 2013). The dashed line is the best linear fit of the 
CH 3 OH maser features of group A (PAchjOH = -49° ± 7°). 


the other maser features we found that - 55°| > \6^ - 55°|. We 
did not detect any circularly polarized maser emission toward 
the region (Py < 1.5%). 
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Fig. 4. View of the CH 3 OH maser features detected around G29.86-0.04 
(Table A.3). Same symbols as in Fig. 2. The assumed velocity of the 
massive SFR region is = -1-101.85 km s“* (de Villiers et al. 2014). 
The two arrows indicate the direction, and not the actual position, of the 
red- and blueshifted lobe of the bipolar outflow ~ - 1 - 6 ° and 

P'^bi™shifted ~ +60°; de Villiers et al. 2014). The dashed lines are the 
best linear fit to the positions of the CFI 3 OH maser features of group A 
and B (PA^HjOH = +8 ±7 and PA®jj^Qpj = -49 ± 5 ). 


Fig. 5. View of the CFI 3 OH maser features detected around 
G35.03-l-0.35 (Table A.4). Same symbols as in Fig. 2. The assumed 
velocity of the massive SFR region is Plsr = -)-51.5 kms * (Paron et 
al. 2012). The dotted line indicates the direction, and not the actual po¬ 
sition, of the bipolar 4.5 yum emission (PA'*-^''™ sj -i-27°; Cyganowski 
et al. 2009) and the two arrows indicate the direction of the red- and 
blueshifted lobes of the *^CO-outflow (Paron et al. 2012). The dashed 
line is the best linear fit of the CH 3 OH maser features (PAchjOH = 
-26° ± 19°). 


4.3. G29.86-0.04 

In Table A.3 and Fig. 4 we report the 18 CH 3 OH maser fea¬ 
tures that we detected in the region. We divided the maser fea¬ 
tures into two groups (A and B), and from a linear fit we find 
that the features of group A are aligned with the redshifted lobe 
of the outflow (PAA^^o^ = + 8 ° ± 7°and ^ + 6 °). 

The five maser features of group B are instead linearly dis¬ 
tributed perpendicularly to the blueshifted lobe of the outflow 
(PA?H 30 h = -49° ± 5°and , +63°), even though 

they are spatially associated with the redshifted lobe (see Fig. B- 
1 of de Villiers et al. 2014). The maser features of group A show 
a velocity gradient, from north (the most blueshifted velocity) 
to south (the most redshifted velocity), the range of which is 
consistent with the velocities of the quiescent emission of '^CO 
(-H96.5kms-'< , < h- 104km s^'; de Villiers et al. 2014). 

^ quiescent ’ ^ 

The velocity range of group B is also consistent with Fqyfg® 
Almost 40% of the CH 3 OH maser features show linearly po¬ 
larized emission (Pi = 1.2% -17%) and only the highest linearly 
polarized feature (i.e., G29.17 in Table A.3) appears to have a 
high saturation degree {T\,AQ. - 6.3 x 10’° K sr). From our 
analysis of the estimated 6 values we determined that the mag¬ 
netic held is perpendicular to all the maser features but G29.14, 
for which 6 - 67°4^i°°. We measured a Zeeman splitting of 
AVz = - 6.6 +1.1 ms“' toward G29.09 {Py = 0.5%; see Fig. 1). 

4.4. G35.03+0.35 

Across a bandwidth that covers a range of velocities be¬ 
tween -6 kms^' and +94 kms^' , we detected 29 6.7- 
GHz CH 3 OH maser features with velocities +41.5 kms“^< 
Fisr <+46.7 kms“*(see Table A.4). No redshifted features were 


detected, as previously reported (Szymczak et al. 2000; Pandian 
et al. 2011). The maser features (Fig. 5) are distributed from 
southeast to northwest (PAchjOh = -26° + 19°) almost perpen¬ 
dicular to the 4.5 yum emission. In Fig. 5 we have drawn the two 
arrows assuming that the bipolar 4.5 /urn emission also traces the 
'^CO large-scale outflow (Cyganowski et al. 2009, Paron et al. 
2012 ). 

Because of the weak 6.7 GHz CH 3 OH maser features, we 
were able to measure linear polarization only toward the bright¬ 
est maser feature G35.19 (Pi = 0.9%), which appears to be un¬ 
saturated. The corresponding 0 angle was +90° + 24°. No circu¬ 
larly polarized emission was detected (Py < 0.8%). 

4.5. G37.43+1.51 

We detected two groups of CH 3 OH maser features, named 
group A and group B in Fig. 6 , separated by 300 mas (~550 au). 
Group A is composed of 14 maser features distributed linearly 
with PAchjOh = -64° + 5° with no clear velocity gradient, as 
already reported by Fujisawa et al. (2014). The velocities of 
group A are consistent with the velocity range of the blueshifted 
lobe of the C'^O-outflow (Lopez-Sepulcre et al. 2010). Maser 
features of group B were not detected before. This group, which 
is located southeast w.r.t. group A, show a velocity range of be¬ 
tween 46 km s“^ and 52 km s“'. Furthermore, an isolated maser 
feature (G37.07; see Table A.5) that cannot be associated with 
either of the two groups is located at 400 mas (~750 au) south 
and 300 mas (~560 au) west from groups A and B, respectively. 

We detected linearly polarized emission from three 
CH 3 OH maser features of group A (Pi = 0.7% - 1.5%), all of 
which have an estimated PbAQ lower than the saturation thresh¬ 
old. The FRTM code estimated that the magnetic held is paral¬ 
lel to all the linear polarization vectors of these features, indeed 
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Fig. 6. View of the CH3OH maser features detected around 
G37.43+1.51 (Table A.5). Same symbols as in Fig. 2. The assumed 
velocity of the massive SFR region is = +44.1 kms“* (Lopez- 

Sepulcre et al. 2010). The two arrows indicate the direction, and not the 
actual position, of the red- and blueshifted lobes of the C'^O-outflow 
Lopez-Sepulcre et al. 2010). The dashed line is the 
best linear fit of the CFlsOFI maser features (PACH 3 OH = -64° ± 5°). 



Fig. 7. View of the CFI 3 OFI maser features detected around the mil¬ 
limeter core MMl of G174.20-0.08 (Table A. 6 ). Same symbols as in 
Fig. 2. The assumed velocity of the massive SFR region is = 

-1.0 kms^' (Zhang et al. 2007). The two arrows indicate the direc¬ 
tion, and not the actual position, of the red- and blueshifted lobes of the 
‘^CO-outflow (PA“2° „ = -40°; Goddi et al. 2011). 

10+ - 55°| < \6- - 55°| (see Col. 14 of Table A.5). No circular 
polarization was measured (Py < 0 . 2 %). 

4.6. G174.20-0.08 

The 14 6.7-GHz CH 3 OH maser features detected toward 
AFGL5142 are shown in Fig. 7. No CH 3 OH maser emission 


Fig. 8. View of the CFI3OFI maser features detected toward Star A in 
the massive SFR G213.70-12.6 (Mon R2-IRS3; Table A.7). Same sym¬ 
bols as in Fig. 2. The assumed velocity of the massive SFR region is 
“ +10.5 kms^' (Dierickx et al. 2015). The arrows indicate 
the direction, and not the actual position, of the red- and blueshifted 
*'^CO(2-l)-outflow powered by IRS 3 = +53°; Dierickx et 

al. 2015). The dashed line is the best linear fit of the CH 3 OH maser 
features (PAchjOH = +63° ± 2°). 

with a peak flux density > 0.9 Jybeam * was detected. Both 
the maser distribution and the velocity range of the masers agree 
with previous observations (e.g., Goddi et al. 2011). We were 
not able to detect at 5 cr either linearly polarized (Pi < 0 .02%) or 
circularly polarized maser emissions (Py < 0.02%). 

4.7. G213.70-12.6 

We detected 20 CH 3 OH maser features that are linearly dis¬ 
tributed from northeast to southwest with PAchjOh = +63° + 
2°(see Fig. 8 ). Because the most western maser features 
(G213.01-G213.04 in Table A.7) were previously undetected 
(Minier et al. 2000), the linear distribution of the maser features 
is now more extended (545 mas; ~450 au). 

Six of the CH 3 OH maser features showed linearly polarized 
emission (Pi = 3.0% - 5.0%), and according to the estimated 
PbAQ three of them are unsaturated (G213.08, G213.12, and 
G213.13). The FRTM code estimated 0 angles greater than 55°, 
indicating that the magnetic held is perpendicular to all the mea¬ 
sured linear polarization vectors. Furthermore, we detected a cir¬ 
cular polarization of 0.6% toward the brightest CH 3 OH maser 
feature (G213.15), which implies a Zeeman splitting of AVz = 
- 6 . 6 + l.Oms^'. 

5. Discussion 

5.1. Magnetic field orientations 

Linear polarization vectors may undergo a rotation when the 
radiation crosses a medium that is immersed in a magnetic 
held. This phenomenon is known as Faraday rotation. Because 
the polarized maser emission may be affected by two of these 
Faraday rotations, the internal ((Fi) and the foreground Faraday 
rotation (<Ff), we briefly determined whether their effects are 
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Table 2. Comparison between the position angle of the magnetic field, CH 3 OH maser distribution, outflows, and linear polarization angles. 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

Source 

Of« 

(x? 

(Kb)'’ 

PAoutflOW 

PACH3OH 

P" 

|PAoutllow ~ (^b)I 

IPACH3OH - Cr>l 

IPACH3OH - PAoutflowl 

ref.* 


(°) 

n 

C) 

(°) 

(°) 


(°) 

(°) 

(°) 


IRAS 20126+4104 

4 

-70± 16 

+20 ± 16 

-65 ± 5' 

+87 ±4 

+0.12 

85± 17 

23 ± 17' 

28 ±6* 

(1), (2) 

G24.78+0.08-A2 

17 

-53 ±2 

+37 ± 2^ 

-40± 15* 

-26± 19 

-0.77 

77 ± 15 

79± 19 

66 ±24 

(3), (4) 

G25.65+1.05 

7 

-80 ±8 

-23±5F 

-15 ±15* 

-49 ± r 

-0.87 

8 ±53 

31 ± 11 

64 ± 17 

(3), (5) 

G29.86-0.04 

17 

+46 ±41 

+82 ± 56^ 

+6 ± 15‘* 

+8 ±7' 

+0.73 

76 ±58 

38 ±42 

14± 17 

(3), (6) 

G35.03+0.35 

8 

-64 ±5 

+26 ± 5^ 

+27 ± 15*-* 

-26 ± 19 

-0.77 

1 ± 16 

38 ±20 

53 ±24 

(3), (7), (8) 

G37.43+1.51 

4 

+90 ±3 

+90 ± 3^ 

-4 ± 15* 

-64 ± 5' 

-0.87 

86 ± 15* 

26 ±6* 

60± 16 

(3), (9) 

G174.20-0.08 

4 

- 

- 

-40 ±15* 

-63 ± 16 

-0.45 

- 

- 

23 ±22 

(3), (10) 

G213.70-12.6-IRS3 

2 

+20 ±5 

-70 ± 5^ 

+53 ±15* 

+63 ± 2 

+0.95 

57± 16* 

43 ±5 

10± 15 

(3), (11) 

From Paper IF . 

Cepheus A 

2 

-57 ± 28 

+30± 19 

+40 ±4 

-79 ±9 

-0.34 

10± 19 

22 ±29 

61 ± 10 

(12) 

W75N-group A 

3 

-13±9 

+77 ±9 

+66 ± 15 

+43 ± 10 

+0.96 

11 ± 18 

56± 14 

23± 18 

(12) 

NGC7538-IRS1 

6 

-30 ± 69 

+67 ± 70 

-40± 10 

+84 ±7 

+0.15 

73 ±71 

66 ±69 

56± 12 

(12) 

W3(OH)-group II 

4 

+21 ±45 

-47 ±44 

- 

-59 ±6 

-0.84 


80 ±45 

- 

(12) 

W51-e2 

12 

+33± 16 

-60 ±21 

-50 ± 20 

+57 ±8 

+0.70 

10 ±29 

24± 18 

73 ±22 

(12) 

IRAS18556+0138 

5 

-2± 11 

+88± 11 

+58 ± 23 

-40 ±2 

-0.99 

30 ±26 

42± 11 

82 ±23 

(12) 

W48 

7 

+23 ±7 

-67 ±7 

- 

+55 ± 10 

+0.70 

- 

78 ± 12 

- 

(12) 

IRAS06058+2138-NIRS1 

4 

+49 ± 47 

-49 ± 52 

-50± 15 

+78 ±7 

+0.64 

1 ±54 

29 ±48 

52± 17 

(12) 

IRAS22272+6358A 

2 

-80± 15 

+9± 15 

-40 ± 15 

-35± 11 

-0.87 

49 ±21 

45± 19 

5± 19 

(12) 

S255-IR 

4 

+36± 12 

-54± 12 

+75 ± 15 

-63 ± 49 

-0.11 

51 ± 19 

81 ±51 

42±51 

(12) 

S231 

4 

+28 ± 49 

-62 ± 49 

-47 ±5 

+28 ±8 

+0.97 

15 ±49 

0±50 

75 ±9 

(12) 

G291.27-0.70 

7 

-32 ±5 

+52 ±5 

- 

-77 ± 14 

- 

- 

45± 15 

- 

(12) 

G305.21+0.21 

9 

-51 ± 14 

28± 14 

- 

+48 ± 23 

- 

- 

81 ±27 

- 

(12) 

G309.92+0.47 

12 

+2 ±56 

-75 ± 56 

- 

+35 ±5 

- 

- 

33 ±56 

- 

(12) 

G316.64-0.08 

3 

-67 ± 36 

+21 ±36 

- 

+34 ± 29 

- 

- 

79 ±46 

- 

(12) 

G335.79+0.17 

8 

+44 ±28 

-41 ±28 

- 

-69 ± 25 

- 

- 

67 ±38 

- 

(12) 

G339.88-1.26 

7 

+77 ± 24 

-12 ±24 

- 

-60 ± 17 

- 

- 

43 ±29 

- 

(12) 

G345.01 + 1.79 

5 

+5 ±39 

-86 ± 39 

- 

+74 ±4 

- 

- 

69 ±39 

- 

(12) 

NGC6334F (central) 

5 

+77 ± 20 

-13 ±20 

+30± 15* 

-41 ± 16 

- 

43 ±25 

62 ±26 

71 ±41 

(12); (13) 

NGC6334F (NW) 

5 

-71 ±20 

+ 19 ±20 

+30± 15* 

-80 ±38 

- 

11 ±25 

9 ±43 

70 ±41* 

(12); (13) 


Notes. Foreground Faraday rotation estimated by using Eq. 3 of Paper I. Because of the high uncertainties of the estimated Of, the angles 
are not corrected for Of. Pearson product-moment correlation coefficient -1 <p< +\\ p = H-l (p = - l)is total positive (negative) correlation, 
p = 0 is no correlation. (1) Surcis et al. (2014); (2) Moscadelli et al. (2011); (3) this work; (4) Beltran et al. (2011); (5) Sanchez-Monge et al. 
(2013); (6) de Villiers et al. (2014); (7) Cyganowski et al. (2009); (8) Paron et al. (2012); (9) Lopez-Sepulcre et al. (2010); (10) Goddi et al. (2011); 
(11) Dierickx et al. (2015); (12) Paper II and references therein; (13) Zhang et al. (2014). We overestimate the errors by considering half of the 
opening angle of the outflow. The differences between the angles are evaluated taking into account that PA = PA ± 180°, ix) = ix) ± 180°, and 
(Ob) = (®b) ± 180°. Before averaging, we use the criterion described in Sect. 3 to estimate the orientation of the magnetic field w.r.t the linear 
polarization vectors. We consider an arbitrary conservative error of 15°. We consider only group A. ® We consider the PA of the redshifted 
lobe of the CO outflow; see Sect. 5.1. ® We assumed PAoutHow = see Sect. 4.4. Here we omit all the notes that are already included in 

Table 2 of Paper II. 


negligible or not. The former, that is, (Pi, can be considered 
negligible as explained in Papers I and II, while Of needs to 
be estimated numerically by using Eq. 3 of Paper I. We find 
that Of ranges between about 2° and 17°, for four sources it is 
within the errors of the measured linear polarization angles (see 
Tables A.1-A.7), and for three sources it is larger. However, Of 
is very uncertain because the errors of some parameters used to 
calculate it cannot be estimated. Therefore we did not correct 
either the (x) angles or the (Ob) angles, but we list Of in Col.2 
of Table 2 for reader judgment. 

We now discuss separately the orientation of the magnetic 
field in the massive SFRs toward which we detected linearly 
polarized CH 3 OH maser emission. 


G24.78+0.08. Taking into account that for G24.23 and G24.52 
the magnetic field is derived to be parallel to the linear polariza¬ 
tion vector, the error-weighted orientation of the magnetic field 
aroundAl andA2is (O^') = h- 39°+42° and(0^2^ = -h 37° + 2°, 
respectively. Although for G24.784-0.08 Of = 17°, the mag¬ 
netic fields in both cores are oriented preferentially along the 
velocity gradient of the toroidal structures (PAai = 4-50° and 
PAa 2 = 4-40°; Beltran et al. 2011), and not along the CO- 
outflow = -40°; Beltran et al. 2011), indicating that 

the magnetic field is possibly located on their surfaces (see 
Fig. 2). Furthermore, in Al the Zeeman-splitting measurements 


are spatially distributed with the negative measurement in the 
northern maser group and the positive measurement in the 
southern maser group (see left panel of Fig. 2). We recall that 
if AVz > 0 , the magnetic field points away from the observer, 
and if AVz < 0 , toward the observer, the magnetic field around 
Al shows a counterclockwise direction that is opposite to the 
rotation of the toroidal structure. This is similar to what Surcis 
et al. (2011a) measured in NGC7538. We also note that from 
our measurements the magnetic field seems to wrap the gas 
along the preferential southeast-northwest direction of star 
formation (Beltran et al. 2011). We measured AVz < 0 toward 
A2, but in this case, because we have only one measurement, 
we cannot determine if the magnetic field behaves similarly to 
the field associated with Al. Unfortunately, we cannot discern 
if the magnetic field is associated directly with the two toroidal 
structures or with the gas that surrounds all the cores. 

G25.65+1.05. We measured an error-weighted orientation of the 
magnetic field of (<1>b) = -23° + 51° by taking into account that 
for G25.02 and G25.06 |0+ - 55°| < |r - 55°| (see Sect. 4.2). 
Therefore, the magnetic field is oriented along the SiO outflow 
(P^omflow = -15°; Sanchez-Monge et al. 2013). For G25.02 
\0^ - 55°| is smaller than \6^ - 55°| by only 1°. This indicates that 
the probability that the magnetic field is parallel to the linear 
polarization vector is not as high as to completely exclude the 
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Fig. 9 . Left: The probability distribution function (PDF, top panel) and the cumulative distribution function (CDF, bottom panel) of the projected 
angle between the PA of the CH3OFI maser distribution and the linear polarization angles (IPAchsOH - Ct)l)- Right: The PDF and the CDF of the 
projected angle between the PA of the CH3OFI maser distribution and the outflow axes (IPAchjOH “ PAoutflowi)- In both panels the dashed line is the 
CDF for random orientation of outflows and magnetic fields, i.e., all angular differences are equally likely. The results of the Kolmogorov-Smirnov 
test are listed in Table 3. 


opposite. However, even if we consider that the magnetic field 
is perpendicular to the linear polarization vector of G25.02, 
we still have that the magnetic field ((Og) = +1° + 37°) is 
preferentially oriented along the outflow. If we now compare our 
measurements, both (Ob) and (Og), with the measurement of 
the magnetic field at arcsecond scale = + 8 ° + 16°; Val- 

lee & Bastien 2000), we find a good agreement within the errors. 

G29.86-0.04. Although Of = 17° is estimated to be large, the 
magnetic field is oriented almost east-west on the plane of the 
sky ((Ob) = -1-82° + 56°), which is consistent with the orienta¬ 
tion of the bent blueshifted lobe of the CO-outflow. However, 
because the magnetic field orientation is estimated from the 
polarized emission of masers that are spatially associated with 
the redshifted lobe of the CO-outflow, we must only compare 
it with the orientation of the redshifted lobe. Therefore, the 
magnetic field is almost perpendicular to it, suggesting that 
perhaps the CH 3 OH masers probe a magnetic field that might be 
twisted around the axis of the redshifted lobe of the CO-outflow. 
However, because we have only one Zeeman-splitting measure¬ 
ment, which indicates that the magnetic field is pointing toward 
the observer, our interpretation is merely speculative. 

G35.03+0.35. We were able to determine the magnetic field 
orientation on the plane of the sky from one linear polarization 
measurement (Ob = h- 26° + 5°). The magnetic field is oriented 
along the 4.5 pm emission and the projection on the plane of the 
sky of the CO-outflow. 

G37.43+1.51. The magnetic field is assumed to be parallel to all 
the linear polarization vectors of the 6.7 GHz CH 3 OH masers 
measured toward group A of this massive SFR, that is, 
(Ob) = 4-90° + 3°. The magnetic field is thus perpendicular 
to the orientation on the plane of the sky of the C**0-outflow 
Lopez-Sepulcre et al. 2010). 


G213.70-12.6. The magnetic field at mas resolution shows an 
orientation on the plane of the sky of (Ob) = -70° + 5°. The 
magnetic field thus appears to be aligned with the linear polar¬ 
ization vector of Star A (PApJ" ^ = - 68 ° + 2°; Simpson et al. 
2013), which the CH 3 OH masers are associated with, but it is 
rotated by about 90° with respect to the magnetic field inferred^ 
from the polarimetric measurements at scales larger than 0".2 
(Yao et al. 1997; Simpson et al. 2013). The magnetic field probed 
by the masers is almost aligned with the large-scale CO-outflow 
detected toward IRS 6 -45°; Xu et al. 2006), but 

it is almost perpendicular to the small-scale '^CO(2-l) outflow 
associated with IRS 3 (PA^^jj^^^ ss 4-53°; Dierickx et al. 2015). 
Similarly to G29.86-0.04, we here speculate that the magnetic 
field in G213.70-12.6 might be twisted along the outflow axis. 

The Zeeman splitting is measured from the circularly polar¬ 
ized spectra of the brightest maser G213.15 (Fig. 1). Because 
G213.15 is assumed to be partially saturated (see Sect. 4.7), the 
circular polarization might be influenced by a non-Zeeman effect 
due to the saturation state of the maser, that is, the rotation of the 
axis of symmetry for the molecular quantum states (e.g., Vlem- 
mings 2008). Although this effect is difficult to quantify, we are 
quite confident that the contribution to Py of this non-Zeeman 
effect is not high enough to invert the S -shape of the V spec¬ 
tra. Otherwise, we would have measured a much higher value of 
TbAQ for G213.15, which is slightly above the saturation thresh¬ 
old of log(TbAQ) = 9.4 log(K sr). Consequently, from the sign 
of the Zeeman splitting, we can conclude that the magnetic field 
is pointing toward the observer. 

5.2. Updated statistical results 

At the midpoint of our project to determine if there exists 
any relation between the morphology of the magnetic field on 
a scale of tens of astronomical unit and the ejecting direc- 

^ If the magnetic field is considered to be perpendicular to the linear 
polarization vectors of the infrared emissions. 
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tion of molecular outflow from massive YSOs, we must up¬ 
date our first statistical results reported in Paper II by adding 
the new magnetic field measurements made around the sources 
discussed in Sect. 5.1 and around IRAS 20126-1-4104 (Surcis 
et al. 2014). Moreover, we also added two of the southern 
sources observed by Dodson & Moriarty (2012) to our analysis: 
NGC6334(central) and NGC6334(NW), which were recently as¬ 
sociated with the blueshifted lobe of a CO-outflow (Zhang et 
al. 2014). Therefore we analyzed the probability distribution 
function (PDF) and the cumulative distribution function (CDF) 
of the projected angles |PAoutflow - <<E>b>I, IPAchjOh - (t>L and 
IPAchjOh - PAoutfiowl; where PAoutfiow is the orientation of the 
large-scale molecular outflow on the plane of the sky, (Ob) is the 
error -weighted orientation of the magnetic field on the plane of 
the sky, PAchjOh is the orientation of the CH 3 OH maser distribu¬ 
tion, and ix) is the error-weighted value of the linear polarization 
angles. Note that although Surcis et al. (2014) determined the 
morphology of the magnetic field around IRAS 20126-1-4104 by 
observing the polarized emission of both 6.7 GHz CH 3 OH and 
22 GHz H 2 O masers, we consider here only the orientation of 
the magnetic field estimated from the CH 3 OH masers. 

We list all the sources of the updated magnetic field total 
sample in Table 2; note that all the angles are the projection on 
the plane of the sky. For the statistical analysis we require the un¬ 
certainties of all the angles. While the errors of PAchsOh and ix) 
are easily determined, the uncertainties of PAoutfiow are unknown 
for all the new sources but IRAS 20126-1-4104. Therefore, as al¬ 
ready done in Paper II, we considered a conservative uncertainty 
of +15°. The uncertainties in Cols. 8 to 10 of Table 2 are equal 

to c+x-y = where x and y are the two angles taken in 

consideration in each column. 

In Figs. 9 and 10 we show the PDF and the CDF of 
IPAchjOh - (t)L IPAchjOh - PAoutfiowl, and IPAoutfiow - (^b)!- The 
results of the Kolmogorov-Smirnov (K-S) test are reported in Ta¬ 
ble 3. We note that the probability that the angles IPAchjOh - ix)\ 
are drawn from a random distribution is now ~80%, which is 
20% higher than what we computed in Paper II. On the other 
hand, the probability for the angles IPAchsOh - PAoutfiowl de¬ 
creases to 34%, which was 60% in Paper II. Note that if more 
than one maser group is detected toward an SFR region, we con¬ 
sider in our analysis the maser group that shows the longest lin¬ 
ear distribution and that is clearly associated with the outflow. 
Even for scattered maser distribution (e.g., G174.20-0.08 MMl) 
we perform a linear fit. 

Although the number of sources for which molecular out¬ 
flows have been detected and for which the orientation of the 
magnetic field has been determined is now twice that of Paper II 
(18 vs. 9), the probability that the distribution of IPAoutffow-^d'B)! 
values are drawn from a random distribution is still 10%. This 
probability confirms our previous conclusion: the magnetic field 
close to the central YSO (10-100 au) is preferentially oriented 
along the outflow axis. 

A more accurate statistical analysis will be presented in the 
last paper of the series when all the sources are observed and 
analyzed. 

6. Summary 

We observed seven massive star-forming regions at 6.7 GHz 
in full polarization spectral mode with the EVN to detect the 
linearly and circularly polarized emission of CH 3 OH masers. 
We detected linearly polarized emission toward all the sources 
but G174.20-0.08 (AEGL5142) and circularly polarized emis- 



Fig. 10. Probability distribution function (PDF, top panel) and the cu¬ 
mulative distribution function (CDF, bottom panel) of the projected an¬ 
gle between the magnetic field and the outflow axes (|PAo„tflow “ (®b)I)- 
The dashed line is the CDF for random orientation of outflows and mag¬ 
netic fields, i.e., all angular differences are equally likely. The results of 
the Kolmogorov-Smimov test are listed in Table 3 

Table 3. Results of the Kolmogorov-Smimov test. 


(1) 

Angle 

(2) 

(3) 

Dh 

(4) 

(5) 

GK-siiP 

IPACH3OH - Ct)l 

21 

0.12 

0.65 

0.79 

IPACH3OH - PAoutfiowl 

19 

0.21 

0.94 

0.34 

IPAoutfiow “ {^b)I 

18 

0.28 

1.22 

0.10 


Notes. N is the number of elements considered in the K.-S. test. 

D is the highest value of the absolute difference between the data 
set, composed of N elements, and the random distribution. 4 is a 
parameter given by 4 = (VY + 0.12 + 0.11/VY) x D. 2k-s(4) = 
2 jg significance level of the K-S test. 


sion toward three sources, G24.78+0.08, G29.86-0.04, and 
G213.70-12.6. By analyzing the polarized emission of the 
masers, we were able to estimate the orientation of the magnetic 
field around seven massive YSOs, considering that G24.78+0.08 
hosts two centers of CH 3 OH masers around each of the YSOs 
Al and A2. The magnetic field is oriented along the outflows 
in two YSOs, it is almost perpendicular to the outflows in four 
YSOs, and in one YSOs (G24.78+0.08 Al) a comparison is not 
possible. Moreover, in G24.78+0.08 Al and A2 the magnetic 
field is oriented along the toroidal structures. From the circularly 
polarized emission of the CH 3 OH masers we measured Zeeman 
splitting toward G24.78+0.08 (both Al and A2), G29.86-0.04, 
and G213.70-12.6. 

We added all the magnetic field measurements made toward 
the YSOs presented in this work to the magnetic field total 
sample, which contains all the massive YSOs observed so far 
in full polarization mode at 6.7 GHz anywhere on the sky. 
Similarly to Paper II, we compared the projected angles between 
magnetic fields and outflows. We still find evidence that the 
magnetic field around massive YSOs are preferentially oriented 
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along the molecular outflows. Indeed, the Kolmogorov-Smirnov 
test still shows a probability of 10% that our distribution of 
angles is drawn from a random distribution. 
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Appendix A: Tables 

In Tables A.1-A.7 we list the parameters of all the 
CH 3 OH maser features detected toward the seven massive star¬ 
forming regions observed with the EVN. The tables are orga¬ 
nized as follows. The name of the feature is reported in Col. 1 
and the group to which they belong in Col. 2. The positions, 
Cols. 3 and 4, refer to the maser feature used for self-calibration. 
The peak flux density, the LSR velocity (Visr), and the FWHM 
(Avl) of the total intensity spectra of the maser features are re¬ 
ported in Cols. 5 to 7. The peak flux density, Visr, and Avl are 
obtained using a Gaussian fit. The mean linear polarization frac¬ 
tion (Pi) and the mean linear polarization angles (x) that are mea¬ 
sured across the spectrum are reported in Cols. 8 and 9. The best¬ 
fitting results obtained by using a model based on the radiative 
transfer theory of methanol masers for T - 1 - Ty = Is"' (Vlem- 
mings et al. 2010, Surds et al. 2011b) are reported in Cols. 10 
(the intrinsic thermal linewidth) and 11 (the emerging brightness 
temperature). The errors were determined by analyzing the full 
probability distribution function. The angle between the mag¬ 
netic field and the maser propagation direction (0, Col. 14) is 
determined by using the observed Pi and the fit emerging bright¬ 
ness temperature. The errors for 0 were also determined by an¬ 
alyzing the full probability distribution function. The value of 0 
in bold indicates that \0^ - 55°| < \0^ - 55°|, that is, the magnetic 
field is assumed to be parallel to the linear polarization vector 
(see Sect. 3). The circular polarization fraction (Py) and the Zee- 
man splitting (AVz) are listed in Cols. 12 and 13. The Zeeman 
splitting is determined by fitting the V Stokes spectra by using 
the best-fitting results (AV) and PbAQ). 
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Table A.l. Parameters of the 6.7 GHz CH 3 OH maser features detected in G24.78+0.08. 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

— 

(10) 

(11) 

(12) 

(13) 

(14) 

Maser 

Group 

RA'^ 

Dec^ 

Peak flux 

Vk, 

Akl 

Pi 



AVi 

TbASl 

Pv 

AVz 

9 



offset 

offset 

density 













(mas) 

(mas) 

(Jy/beam) 

(km/s) 

(km/s) 

(%) 

(°) 


(km/s) 

(log K sr) 

(%) 

(m/s) 

0 

G24.01 

A2 

-1542.297 

746.773 

0.092 ± 0.005 

108.57 

0.25 

- 

- 


- 

- 

- 

- 

- 

G24.02 

A2 

-1449.196 

553.302 

0.222 ± 0.004 

106.59 

0.28 

- 

- 


- 

- 

- 

- 

- 

G24.03 

A2 

-1416.923 

510.347 

0.031 ±0.004 

105.41 

0.20 

- 

- 


- 

- 

- 

- 

- 

G24.04 

A2 

-1410.906 

1117.249 

0.115±0.014 

113.18 

0.47 

- 

- 


- 

- 

- 

- 

- 

G24.05 

A2 

-1410.054 

1112.544 

0.108 ±0.006 

113.00 

0.34 

- 

- 


- 

- 

- 

- 

- 

G24.06 

A2 

-1404.405 

1082.697 

0.393 ± 0.006 

113.00 

0.25 

- 

- 


- 

- 

- 

- 

- 

G24.07 

A2 

-1396.614 

1063.488 

0.697 ±0.019 

111.95 

0.24 

- 

- 


- 

- 

- 

- 

- 

G24.08 

A2 

-1395.337 

1044.547 

0.086 ±0.010 

113.13 

0.35 

- 

- 


- 

- 

- 

- 

- 

G24.09 

A2 

-1380.307 

1059.353 

0.437 ±0.014 

112.43 

0.30 

- 

- 


- 

- 

- 

- 

- 

G24.10 

A2 

-1391.533 

1035.954 

4.466 ±0.019 

111.95 

0.27 

1.3±0.1 

-51 ± 

3 

^•^- 0.2 

8.7(2? 

- 

- 

83(2, 

G24.il 

A2 

-1387.531 

1015.326 

0.233 ± 0.006 

111.64 

0.21 

- 

- 


- 

- 

- 

- 

- 

G24.12 

A2 

-1383.656 

1059.686 

0.524 ±0.018 

111.95 

0.70 

- 

- 


- 

- 

- 

- 

- 

G24.13 

A2 

-1358.919 

984.020 

0.261 ±0.006 

111.64 

0.28 

- 

- 


- 

- 

- 

- 

- 

G24.14 

A2 

-1341.548 

1037.778 

0.118 ±0.004 

110.89 

1.68 

- 

- 


- 

- 

- 

- 

- 

G24.15 

A2 

-1332.025 

1310.781 

0.101 ±0.006 

110.10 

0.65 

- 

- 


- 

- 

- 

- 

- 

G24.16 

A2 

-1324.574 

694.642 

13.699 ± 0.036 

110.41 

0.32 

1.1±0.1 

-55 ± 

5 

1.7fS;2 

8.6(2? 

0.3 

-4.0 ±0.8 


G24.17 

A2 

-1310.680 

726.521 

0.234 ± 0.006 

109.58 

0.27 

- 

- 


- 

- 

- 

- 

- 

G24.18 

A2 

-1307.118 

705.516 

0.044 ± 0.005 

109.84 

0.13 

- 

- 


- 

- 

- 

- 

- 

G24.19 

A2 

-1216.628 

836.181 

0.054 ± 0.004 

110.98 

0.22 

- 

- 


- 

- 

- 

- 

- 

G24.20 

A2 

-1155.445 

997.989 

2.826 ±0.019 

112.12 

0.37 

1.0 ±0.1 

-53 ± 

8 

1.8t^ 

8.6(2i 

- 

- 

82(2, 

G24.21 

A1 

-536.18 

-377.044 

0.088 ± 0.004 

108.65 

0.27 

- 

- 


- 

- 

- 

- 

- 

G24.22 

A1 

-482.635 

-407.791 

0.255 ± 0.004 

107.21 

0.21 

- 

- 


- 

- 

- 

- 

- 

G24.23 

A1 

-442.953 

-405.630 

9.791 ±0.025 

107.69 

0.35 

1.3 ±0.2 

-74 ± 

2 

1 9+0.2 

0 O+0.3 
°-°-0.2 

0.5 

+7.8 ± 1.4 


G24.24 

A1 

-438.171 

-331.857 

1.029 ±0.025 

107.78 

0.30 

- 

- 


- 

- 

- 

- 

- 

G24.25 

A1 

-438.157 

-330.956 

0.497 ±0.019 

107.86 

1.90 

- 

- 


- 

- 

- 

- 

- 

G24.26 

G24.27 

A1 

A1 

-433.147 

-426.504 

-406.664 

-404.887 

5.128 ±0.025 
1.564 ±0.009 

107.78 

108.30 

0.29 

0.31 

1.9 ±0.3 
1.3±0.1 

-79 ± 
-64 ± 

2 

4 

1 4+0.2 
-Q.3 

1 5+^-2 
-0.3 

8.9(2i 

8.8(2? 

— 

— 

83+ 

83+ 

G24.28 

A1 

-408.963 

-408.211 

2.670 ±0.018 

108.00 

0.24 

2.2 ±0.4 

-76 ± 

2 

1 

-0.3 

9.0(2? 

- 

- 

81+ 

G24.29 

A1 

-399.852 

-404.841 

1.058 ±0.011 

108.22 

0.25 

3.5 ± 0.6 

-69 ± 

3 


9 T+0.3 

- 

- 

87(?, 

G24.30 

A1 

-382.197 

-403.043 

2.013 ±0.017 

108.00 

0.26 

- 

- 


- 

- 

- 

- 

- 

G24.31 

A1 

-381.913 

-402.766 

2.081 ±0.016 

108.08 

0.27 

0.8±0.1 

-87 ± 

8 

1 4 + 0.2 

8.5(»? 

- 

- 

84+6 

G24.32 

A1 

-375.753 

-401.033 

2.644 ±0.017 

108.04 

0.33 

2.2 ±0.4 

-74 ± 

3 

1 1^^^ 
'- 0.3 

9.0^4 

- 

- 

8i£ 

G24.33 

A1 

-317.792 

-414.340 

0.052 ± 0.005 

107.07 

0.29 

- 

- 


- 

- 

- 

- 

- 

G24.34 

A1 

-287.293 

-413.178 

0.059 ± 0.004 

106.90 

0.21 

- 

- 


- 

- 

- 

- 

- 

G24.35 

A1 

-255.928 

115.137 

0.092 ± 0.004 

116.91 

0.54 

- 

- 


- 

- 

- 

- 

- 

G24.36 

A1 

-177.289 

-116.315 

0.049 ± 0.004 

108.83 

0.23 

- 

- 


- 

- 

- 

- 

- 

G24.37 

A1 

-150.438 

101.504 

0.092 ± 0.004 

115.06 

0.30 

- 

- 


- 

- 

- 

- 

- 

G24.38 

A1 

-142.022 

78.130 

2.812 ±0.006 

114.80 

0.41 

1.8 ±0.5 

-70 ± 

2 

- 

- 

- 

- 

- 

G24.39 

A1 

-128.156 

-76.137 

0.500 ±0.019 

110.28 

0.64 

- 

- 


- 

- 

- 

- 

- 

G24.40 

A1 

-108.358 

107.273 

0.306 ± 0.022 

113.26 

0.32 

- 

- 


- 

- 

- 

- 

- 

G24.41 

A1 

-85.948 

98.951 

1.575 ±0.010 

114.05 

0.31 

- 

- 


- 

- 

- 

- 

- 

G24.42 

A1 

-80.442 

173.600 

0.288 ±0.016 

114.19 

0.30 

- 

- 


- 

- 

- 

- 

- 

G24.43 

A1 

-76.652 

99.772 

11.483 ±0.016 

114.23 

0.33 

1.9 ±0.1 

-71 ± 

2 

1 6+^-2 
^•"- 0.3 

00 

1 + 

D p 

0.7 

- 9.7 ± 1.6 

88!?. 

G24.44 

A1 

-58.415 

-26.639 

0.101 ±0.006 

113.70 

0.43 

- 

- 


- 

- 

- 

- 

- 

G24.45 

A1 

-54.370 

-26.016 

0.175 ±0.008 

113.09 

0.95 

- 

- 


- 

- 

- 

- 

- 

G24.46 

A1 

-37.638 

-28.224 

0.697 ±0.012 

114.32 

0.49 

- 

- 


- 

- 

- 

- 

- 

G24.47 

A1 

-36.900 

-27.994 

0.739 ± 0.006 

113.70 

0.40 

- 

- 


- 

- 

- 

- 

- 

G24.48 

A1 

-32.259 

-24.931 

0.143 ±0.004 

115.06 

0.25 

- 

- 


- 

- 

- 

- 

- 

G24.49 

A1 

-28.796 

108.014 

0.136 ±0.006 

114.71 

0.27 

- 

- 


- 

- 

- 

- 

- 

G24.50 

A1 

-24.822 

-36.330 

0.541 ±0.008 

113.09 

0.28 

- 

- 


- 

- 

- 

- 

- 

G24.51 

A1 

-17.499 

-56.006 

0.829 ± 0.008 

114.41 

0.31 

- 

- 


- 

- 

- 

- 

- 

G24.52 

A1 

0 

0 

32.970 ± 0.032 

113.40 

1.72 

1.5 ±0.4 

-71 ± 

2 

1 c+ 0.2 

0 O+0.5 
-0.3 

0.3 

-4.6 ±0.7 


G24.53 

A1 

59.167 

24.964 

0.047 ± 0.005 

112.78 

0.25 

- 

- 


- 

- 

- 

- 

- 


Notes. The reference position is a 20 oo = +18’'36“12*.563 and ^2000 = -07°12'10"787 (see Sect. 3). 
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Table A.2. Parameters of the 6.7 GHz CH 3 OH maser features detected in G25.65+1.05. 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 


(10) 

(11) 

(12) 

(13) 

(14) 

Maser 

Group 

RA'’ 

Dec“ 

Peak flux 


Avl 

A 

A 


AVi 

7-bAn 

Py 

AVz 

9 



offset 

offset 

density 













(mas) 

(mas) 

(Jy/beam) 

(km/s) 

(km/s) 

(%) 

0 


(km/s) 

(log K sr) 

(%) 

(m/s) 

0 

G25.01 

A 

-19.690 

15.911 

0.338 ±0.014 

41.56 

0.11 

- 

- 


- 

- 

- 

- 

- 

G25.02 

A 

0 

0 

6.924 ±0.124 

41.78 

0.31 

0.4 ±0.2 

-76 ± 

11 

1 6+^‘ 
- 0.2 

8-27:6 

- 

- 

72-11 

G25.03 

A 

6.886 

-0.382 

3.860 ±0.126 

42.00 

0.17 

- 

- 


- 

- 

- 

- 

- 

G25.04 

A 

7.170 

-0.828 

4.023 ±0.108 

42.04 

0.17 

- 

- 


- 

- 

- 

- 

- 

G25.05 

A 

16.873 

-46.883 

0.311 ±0.012 

41.52 

0.11 

- 

- 


- 

- 

- 

- 

- 

G25.06 

A 

18.324 

-13.035 

24.318 ±0.142 

41.82 

0.21 

0.3 ±0.2 

86 ±29 


8-i!l:i 

- 

- 


G25.07 

A 

25.865 

-21.118 

0.045 ± 0.004 

42.66 

0.21 

- 

- 


- 

- 

- 

- 

- 

G25.08 

A 

34.771 

-26.409 

0.448 ±0.012 

41.52 

0.12 

- 

- 


- 

- 

- 

- 

- 

G25.09 

A 

37.787 

-23.865 

0.127 ±0.003 

42.57 

0.22 

- 

- 


- 

- 

- 

- 

- 

G25.10 

A 

42.823 

-33.508 

0.213 ±0.004 

42.66 

0.19 

- 

- 


- 

- 

- 

- 

- 

G25.ll 

A 

43.876 

-64.026 

2.455 ±0.108 

42.04 

0.23 

1.3 ±0.8 

-87 ± 

21 

< 0.5 

8.7-^“-2 

- 

- 

90+g 

G25.12 

A 

51.103 

-23.075 

0.126 ±0.003 

42.75 

0.24 

- 

- 


- 

- 

- 

- 

- 

G25.13 

A 

55.030 

-42.839 

6.113±0.124 

41.78 

0.19 

0.8 ±0.3 

-73 ± 

10 

o.9f!!? 

8.5+;-* 

- 

- 

77 - 1-13 

G25.14 

A 

62.143 

-44.819 

3.748 ±0.108 

42.04 

0.20 

0.9 ±0.2 

-82 ± 

15 


8-5!J:i 

- 

- 

75+15 

G25.15 

A 

63.480 

-41.443 

0.031 ±0.002 

42.97 

0.21 

- 

- 


- 

- 

- 

- 

- 

G25.16 

A 

67.236 

-76.904 

0.016 ± 0.002 

43.76 

0.21 

- 

- 


- 

- 

- 

- 

- 

G25.17 

A 

80.809 

-62.267 

0.041 ±0.003 

42.70 

0.17 

- 

- 


- 

- 

- 

- 

- 

G25.18 

A 

84.252 

-90.500 

0.043 ± 0.002 

43.80 

0.22 

- 

- 


- 

- 

- 

- 

- 

G25.19 

B 

402.507 

-166.119 

0.104 ±0.003 

40.64 

0.27 

- 

- 


- 

- 

- 

- 

- 

G25.20 

B 

423.904 

-168.041 

0.237 ± 0.003 

40.38 

0.29 

- 

- 


- 

- 

- 

- 

- 

G25.21 

B 

428.827 

-92.293 

0.038 ± 0.003 

40.73 

0.23 

- 

- 


- 

- 

- 

- 

- 

G25.22 

B 

442.399 

-214.329 

0.106 ±0.003 

41.03 

0.25 

- 

- 


- 

- 

- 

- 

- 

G25.23 

B 

444.249 

-228.867 

0.094 ± 0.003 

40.99 

0.24 

- 

- 


- 

- 

- 

- 

- 


Notes. The reference position is a 20 oo = +18’'34“20“.900 and S 2000 = -05°59'42"098 (see Sect. 3). 


Table A.3. Parameters of the 6.7 GHz CH 3 OH maser features detected in G29.86-0.04. 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

Maser 

Group 

RA^ 

Dec*^ 

Peak flux 


Avl 

Pi 


AVi 

7-bAn 

Pv 

AVz 

9 



offset 

offset 

density 












(mas) 

(mas) 

(Jy/beam) 

(km/s) 

(km/s) 

(%) 

(°) 

(km/s) 

(log K sr) 

(%) 

(m/s) 

n 

G29.0I 

B 

-188.953 

-153.172 

0.150 ±0.003 

103.68 

0.24 

- 

- 

- 

- 

- 

- 

- 

G29.02 

B 

-183.009 

-160.760 

0.971 ±0.015 

102.05 

0.29 

- 

- 

- 

- 

- 

- 

- 

G29.03 

B 

-157.861 

-185.738 

0.275 ± 0.003 

103.85 

0.32 

- 

- 

- 

- 

- 

- 

- 

G29.04 

B 

-145.858 

-197.453 

2.423 ± 0.004 

104.29 

0.92 

1.24 ±0.03 

67 ±3 

<0.5 

8.9t| 

- 

- 

79-1 

G29.05 

B 

-132.541 

-200.275 

0.053 ± 0.003 

103.06 

0.41 

- 

- 

- 


- 

- 

- 

G29.06 

A 

-11.774 

-50.191 

0.080 ± 0.003 

102.49 

0.24 

- 

- 

- 

- 

- 

- 

- 

G29.07 

A 

-3.715 

-32.043 

3.047 ±0.006 

102.23 

0.20 

3.9 ± 0.5 

19 ±2 

0 . 9-11 

Q a+O-S 

- 

- 

83fJ, 

G29.08 

A 

-3.315 

-20.612 

1.593 ±0.016 

102.01 

0.18 

- 

- 

- 

- 

- 

- 

- 

G29.09 

A 

-3.258 

41.101 

19.262 ±0.018 

100.39 

0.32 

1.6 ±0.3 

11 ±2 



0.5 

-6.6 ± 1.1 

74/!^ 

G29.10 

A 

0 

0 

32.456 ±0.015 

101.88 

0.35 

1.6 ±0.5 

-74 ± 1 

1-8^3 

8.9ii 

- 

- 

73!« 

G29.il 

A 

0.343 

58.858 

0.182 ±0.003 

99.95 

1.69 

- 

- 

- 

- 

- 

- 

- 

G29.12 

A 

0.857 

54.225 

9.863 ±0.009 

100.82 

0.43 

2.1 ±0.3 

25 ±8 

<0.5 

9 O+o.i 

- 

- 

90/1^ 

G29.13 

A 

1.429 

141.063 

0.293 ± 0.003 

99.77 

0.36 

- 

- 

- 


- 

- 


G29.14 

A 

2.286 

7.941 

24.431 ±0.022 

101.48 

0.25 

3.4 ±2.0 

26 ± 1 

1 

^•^- 0.3 

Q 9+0.4 

- 

- 

67-* 

G29.15 

A 

3.372 

145.423 

0.196 ±0.003 

98.89 

0.23 

- 

- 

- 

- 

- 

- 

- 

G29.16* 

A 

5.430 

20.873 

5.179±0.018 

101.26 

0.23 

17.0 ±2.0 

21 ±2 

< 0.5 

10.8f“ 

- 

- 

90!‘ 

G29.17 

A 

13.431 

135.044 

0.083 ± 0.003 

99.03 

3.71 

- 

- 

- 


- 

- 


G29.18 

A 

14.460 

126.231 

0.099 ± 0.003 

99.11 

0.44 

- 

- 

- 

- 

- 

- 

- 


Notes. The reference position is Q '2000 = -i-18''45"'59.572'* and ^2000 = -02°45'01''573 (see Sect. 3). Because of the degree of saturation, 
TbASl is underestimated, AV\ and 6 are overestimated. 
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Table A.4. Parameters of the 6.7 GHz CH 3 OH maser features detected in G35.03+0.35. 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

Maser 

Group 

RA'’ 

Dec" 

Peak flux 


Avl 

A 


AVi 

TbACl 

Pv 

AVz 

9 



offset 

offset 

density 












(mas) 

(mas) 

(Jy/beam) 

(km/s) 

(km/s) 

(%) 

n 

(km/s) 

(log K sr) 

(%) 

(m/s) 

n 

G35.01 


78.111 

64.224 

0.112 ±0.005 

44.88 

0.24 

- 

- 

- 

- 

- 

- 

- 

G35.02 


67.948 

-245.222 

0.061 ±0.006 

46.63 

0.22 

- 

- 

- 

- 

- 

- 

- 

G35.03 


40.077 

74.663 

0.026 ± 0.005 

45.10 

0.23 

- 

- 

- 

- 

- 

- 

- 

G35.04 


37.305 

-208.235 

0.042 ± 0.007 

45.84 

0.31 

- 

- 

- 

- 

- 

- 

- 

G35.05 


31.409 

-140.110 

0.065 ± 0.005 

45.62 

0.32 

- 

- 

- 

- 

- 

- 

- 

G35.06 


25.080 

-138.308 

0.530 ± 0.007 

45.89 

0.38 

- 

- 

- 

- 

- 

- 

- 

G35.07 


18.546 

-114.414 

1.164 ±0.011 

44.26 

0.16 

- 

- 

- 

- 

- 

- 

- 

G35.08 


14.840 

-224.552 

0.055 ± 0.008 

44.62 

0.41 

- 

- 

- 

- 

- 

- 

- 

G35.09 


14.658 

-130.545 

0.035 ± 0.007 

45.93 

0.30 

- 

- 

- 

- 

- 

- 

- 

G35.10 


14.657 

-59.593 

0.087 ± 0.004 

42.99 

0.16 

- 

- 

- 

- 

- 

- 

- 

G35.il 


14.070 

-244.907 

0.099 ± 0.009 

44.09 

0.11 

- 

- 

- 

- 

- 

- 

- 

G35.12 


13.540 

-120.426 

1.062 ±0.013 

44.48 

0.20 

- 

- 

- 

- 

- 

- 

- 

G35.13 


12.015 

-68.813 

0.226 ±0.010 

43.89 

0.46 

- 

- 

- 

- 

- 

- 

- 

G35.14 


11.052 

-169.041 

0.157 ±0.010 

44.22 

0.15 

- 

- 

- 

- 

- 

- 

- 

G35.15 


9.444 

-61.880 

0.182 ±0.012 

44.48 

0.17 

- 

- 

- 

- 

- 

- 

- 

G35.16 


9.048 

-109.791 

0.088 ± 0.009 

44.18 

0.24 

- 

- 

- 

- 

- 

- 

- 

G35.17 


4.780 

-111.321 

0.197 ±0.011 

44.31 

0.20 

- 

- 

- 

- 

- 

- 

- 

G35.18 


1.635 

-182.730 

0.047 ± 0.004 

42.60 

0.25 

- 

- 

- 

- 

- 

- 

- 

G35.19 


0 

0 

2.252 ±0.010 

43.87 

0.46 

0.9 ±0.1 

-64 ±5 

2 

^•^-0.3 

Q c-i-0.3 

- 

- 


G35.20 


-0.201 

-102.259 

0.084 ± 0.005 

43.47 

1.12 

- 

- 

- 

- 

- 

- 

- 

G35.21 


-0.576 

-199.854 

0.087 ± 0.006 

43.34 

0.20 

- 

- 

- 

- 

- 

- 

- 

G35.22 


-1.651 

-173.355 

0.052 ± 0.006 

42.16 

0.24 

- 

- 

- 

- 

- 

- 

- 

G35.23 


-9.665 

-18.393 

0.051 ±0.005 

42.38 

0.47 

- 

- 

- 

- 

- 

- 

- 

G35.24 


-11.233 

-55.342 

0.370 ± 0.005 

42.55 

0.26 

- 

- 

- 

- 

- 

- 

- 

G35.25 


-13.312 

2.804 

0.108 ±0.006 

41.59 

0.28 

- 

- 

- 

- 

- 

- 

- 

G35.26 


-14.343 

-46.606 

0.510 ±0.007 

42.20 

0.21 

- 

- 

- 

- 

- 

- 

- 

G35.27 


-24.309 

-35.076 

0.087 ± 0.004 

42.07 

0.25 

- 

- 

- 

- 

- 

- 

- 

G35.28 


-54.673 

76.618 

0.284 ± 0.006 

46.15 

0.44 

- 

- 

- 

- 

- 

- 

- 

G35.29 


-57.559 

68.207 

0.598 ± 0.006 

46.33 

0.38 

- 

- 

- 

- 

- 

- 

- 


Notes. The reference position is 0:2000 = +18’'54“00“.660 and S 2000 = +02°01'18"551 (see Sect. 3). 
Table A.5. Parameters of the 6.7 GHz CH 3 OH maser features detected in G37.43+1.51. 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

Maser 

Group 

RA" 

Dec" 

Peak flux 

Vkr 

Avl 

Pi 


AVi 

TbAn 

Py 

AVz 

9 



offset 

offset 

density 












(mas) 

(mas) 

(Jy/beam) 

(km/s) 

(km/s) 

(%) 

(°) 

(km/s) 

(log K sr) 

(%) 

(m/s) 

n 

G37.01 

A 

-68.621 

18.314 

0.147 ±0.004 

40.65 

0.21 

- 

- 

- 

- 

- 

- 

- 

G37.02 

A 

-56.601 

15.614 

0.522 ± 0.006 

40.79 

0.20 

- 

- 

- 

- 

- 

- 

- 

G37.03 

A 

-25.050 

-2.205 

0.525 ± 0.004 

42.23 

0.21 

- 

- 

- 

- 

- 

- 

- 

G37.04 

A 

-12.561 

4.368 

17.807 ± 0.064 

41.58 

0.23 

1.5 ±0.4 

89± 18 

^•^- 0.3 

8.8f>:5 

8.7+“ 

- 

- 

71+5 

'^-48 

G37.05 

A 

0 

0 

44.137 ±0.056 

41.36 

0.23 

1.1 ± 1.0 

-87 ± 15 

13 + 0.1 

- 

- 


G37.06 

A 

16.738 

-5.684 

17.809 ± 0.097 

41.09 

0.17 

0.7± 1.1 

87± 14 

0 ^ 
^•^- 0.1 

8.5i| 

- 

- 

G37.07 

- 

22.127 

-412.994 

0.071 ±0.003 

44.60 

0.21 

- 

- 

- 

- 

- 

- 

- 

G37.08 

A 

34.773 

-28.278 

0.133 ±0.003 

40.43 

0.20 

- 

- 

- 

- 

- 

- 

- 

G37.09 

A 

46.621 

-30.720 

0.417 ± 0.004 

40.70 

0.27 

- 

- 

- 

- 

- 

- 

- 

G37.10 

A 

55.475 

-36.919 

0.096 ± 0.004 

42.14 

0.16 

- 

- 

- 

- 

- 

- 

- 

G37.il 

A 

65.854 

-90.763 

0.125 ±0.009 

41.88 

0.27 

- 

- 

- 

- 

- 

- 

- 

G37.12 

A 

86.156 

-50.171 

13.038 ±0.103 

41.14 

0.24 

- 

- 

- 

- 

- 

- 

- 

G37.13 

A 

89.663 

-40.821 

0.379 ± 0.005 

41.97 

0.17 

- 

- 

- 

- 

- 

- 

- 

G37.14 

A 

101.995 

-61.619 

0.597 ±0.019 

40.92 

0.13 

- 

- 

- 

- 

- 

- 

- 

G37.15 

A 

105.916 

-44.788 

3.430 ±0.072 

41.44 

0.20 

- 

- 

- 

- 

- 

- 

- 

G37.16 

B 

315.938 

-401.817 

0.032 ± 0.003 

49.87 

0.35 

- 

- 

- 

- 

- 

- 

- 

G37.17 

B 

326.346 

-381.825 

0.079 ± 0.002 

51.41 

0.42 

- 

- 

- 

- 

- 

- 

- 

G37.18 

B 

336.867 

-376.862 

0.057 ± 0.002 

50.14 

0.29 

- 

- 

- 

- 

- 

- 

- 

G37.19 

B 

371.484 

-331.863 

0.166 ±0.004 

46.05 

0.25 

- 

- 

- 

- 

- 

- 

- 


Notes. The reference position is 0:2000 = +18’'54“ 14^.229 and ^2000 = +04°41'41''138 (see Sect. 3). 
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Table A.6. Parameters of the 6.7 GHz CH3OH maser features detected in G174.20-0.08. 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

Maser 

Group 

RA'’ 

Dec" 

Peak flux 


Akl 

A 


AVi 

GAH 

Pv 

AVz 

9 



offset 

offset 

density 












(mas) 

(mas) 

(Jy/beam) 

(km/s) 

(km/s) 

(%) 

0 

(km/s) 

(log K sr) 

(%) 

(m/s) 

n 

G174.01 


-30.813 

-17.113 

0.047 ± 0.005 

0.43 

0.30 

- 

- 

- 

- 

- 

- 

- 

G 174.02 


-25.677 

-11.940 

0.047 ± 0.007 

1.49 

0.13 

- 

- 

- 

- 

- 

- 

- 

G174.03 


-20.827 

-49.076 

0.057 ± 0.005 

1.88 

0.31 

- 

- 

- 

- 

- 

- 

- 

G 174.04 


-16.548 

-9.541 

0.095 ± 0.008 

1.66 

0.21 

- 

- 

- 

- 

- 

- 

- 

G 174.05 


-12.220 

11.673 

0.043 ± 0.006 

1.53 

0.17 

- 

- 

- 

- 

- 

- 

- 

G 174.06 


-10.271 

-6.065 

0.159 ±0.007 

1.71 

0.17 

- 

- 

- 

- 

- 

- 

- 

G 174.07 


-3.329 

-91.538 

0.105 ±0.005 

3.72 

0.34 

- 

- 

- 

- 

- 

- 

- 

G174.08 


0 

0 

0.853 ±0.006 

1.53 

0.24 

- 

- 

- 

- 

- 

- 

- 

G 174.09 


15.454 

-82.115 

0.029 ± 0.005 

3.72 

0.22 

- 

- 

- 

- 

- 

- 

- 

G174.10 


19.115 

76.355 

0.046 ± 0.009 

5.04 

0.28 

- 

- 

- 

- 

- 

- 

- 

G174.il 


54.778 

-90.572 

0.617 ±0.010 

5.00 

0.29 

- 

- 

- 

- 

- 

- 

- 

G174.12 


73.560 

-80.715 

0.054 ±0.010 

5.00 

0.28 

- 

- 

- 

- 

- 

- 

- 

G174.13 


84.734 

-48.252 

0.135 ±0.005 

2.19 

0.19 

- 

- 

- 

- 

- 

- 

- 

G174.14 


101.234 

-78.148 

0.058 ±0.005 

3.77 

0.20 

- 

- 

- 

- 

- 

- 

- 


Notes. The reference position is 0:2000 = +05’'30“481020 and ^2000 = +33°47'54"611 (see Sect. 3). 


Table A.7. Parameters of the 6.7 GHz CH 3 OH maser features detected in G213.70-12.6 (Mon R2-IRS 3/Star A). 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

Maser 

Group 

RA" 

Dec" 

Peak flux 

Vi» 

Akl 

Pi 


AVi 

TbAn 

Pv 

AVz 

9 



offset 

offset 

density 












(mas) 

(mas) 

(Jy/beam) 

(km/s) 

(km/s) 

(%) 

(°) 

(km/s) 

(log K sr) 

(%) 

(m/s) 

n 

G213.01 


-444.918 

-209.568 

0.111 ±0.006 

10.11 

0.17 

- 

- 

- 

- 

- 

- 

- 

G213.02 


-436.616 

-206.261 

0.098 ± 0.005 

10.07 

0.33 

- 

- 

- 

- 

- 

- 

- 

G213.03 


-334.883 

-188.515 

0.105 ±0.005 

9.98 

0.15 

- 

- 

- 

- 

- 

- 

- 

G213.04 


-304.061 

-187.214 

0.469 ± 0.005 

9.94 

0.18 

- 

- 

- 

- 

- 

- 

- 

G213.05 


-130.791 

-79.498 

3.650 ±0.015 

10.99 

0.21 

3.0 ±0.4 

17 ±2 


A O+0.4 
^•^-0.1 

- 

- 

77+13 

' '-37 

G213.06 


-124.479 

-70.560 

0.570 ±0.015 

10.99 

0.19 

- 

- 

- 

- 

- 

- 

- 

G213.07 


-112.196 

-83.591 

1.139±0.015 

10.99 

0.21 

- 

- 

- 

- 

- 

- 

- 

0213.08* 


-72.276 

-71.392 

49.360 ± 0.06 

10.68 

0.24 

5.0 ± 1.4 

51 ±49 


^•”- 0.5 

- 

- 

71+19 
' ^-27 

G213.09 


-67.784 

-97.511 

1.750 ±0.050 

10.77 

0.18 

- 

- 

- 

- 

- 

- 

- 

G213.10 


-53.624 

-42.149 

0.562 ±0.006 

11.21 

0.21 

- 

- 

- 

- 

- 

- 

- 

G213.il 


-49.189 

-52.086 

0.055 ± 0.004 

12.04 

0.59 

- 

- 

- 

- 

- 

- 

- 

G213.12 


-35.882 

-0.412 

2.470 ± 0.060 

10.68 

0.26 

3.1 ± 0.3 

15±6 

< 0.5 

A O+O.l 

-Q1 

Q 4+0.3 
^ -0.1 

- 

- 

90/13 

G213.13 


-33.721 

-2.464 

5.355 ±0.007 

11.39 

0.26 

4.3 ± 0.3 

23 ±3 

11 + 0.5 
- 0.4 

- 

- 

85/*, 

G213.14 


-11.601 

-12.871 

0.173 ±0.005 

12.31 

0.33 

- 

- 

- 

- 

~ 

- 

- 

G213.15* 


0 

0 

91.580 ±0.050 

12.57 

0.24 

4.6 ±0.4 

16±7 

1 1'-oi 

- 0.4 

9 5+0-^ 

0.6 

-6.6 ± 1.0 

82/?, 

G213.16 


9.269 

-39.421 

0.461 ±0.011 

13.23 

0.25 

- 

- 

- 

- 

- 

- 

- 

G213.17* 


16.377 

13.481 

4.603 ± 0.008 

13.40 

0.25 

3.4 ± 0.3 

21 ±2 

1 1'-o^ 
- 0.3 

9 5+0-^ 

- 

- 

79+18 
' -29 

G213.18 


22.860 

25.299 

0.139 ±0.004 

13.71 

0.27 

- 

- 

- 

- 

- 

- 

- 

G213.19 


28.376 

33.211 

1.702 ±0.009 

13.36 

0.26 

- 

- 

- 

- 

- 

- 

- 

G213.20 


35.086 

49.870 

0.054 ± 0.004 

12.00 

0.67 

- 

- 

- 

- 

- 

- 

- 


Notes. The reference position is Q '2000 = +06’’07“47*.860 and ^2000 = -06°22'56'.'626 (see Sect. 3). Because of the degree of the saturation 
TbAn is underestimated, AV/ and 9 are overestimated. 








